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Enhancement of de Haas-van Alphen Oscillation due to Spin in the 
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The effects of the Zeeman term on the de Haas- van Alphen osciUation is studied in the mag- 
netic breakdown system. We find that the amplitude of the oscillation with the frequencies of 
fp+fa and fp+2fa are enhanced by the Zeeman term, while they are expected to be reduced 
in the semiclassical theory. A possible interpretation of the experiments in organic conductors 
is discussed. 
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§1. Introduction 

The magnetic breakdownEI is observed in the 
Shubnikov-de Haas oscillation (for the transport quan- 
tity) and the de Haas-van Alphen oscillation (for the 
thermodynamics quantity). These oscillations are so- 
called the quantum magnetic oscillation, which is at- 
tributed to the Landau quantization in the closed 
orbit. In quasi- two-dimensional organic conductors 
such as K-(BEDT-TTF)2Cu(NCS)2 and a-(BEDT- 
TTF)2MHg(SCN)4 (M=K, Rb, Tl and NH4), many ex- 
periments for the magnetic breakdown are done, since 
those materials have the suitable Fermi surface (Fig. 1) 
for the study of the magnetic breakdownBtr In this case, 
the electron orbital motion due to magnetic fields {H) 
applied perpendicular to kx-ky plane is confined to the 
small closed orbit named as the a orbit and a pair of open 
orbits at low fields as shown in Fig. 1. Then, this closed 
orbital motion leads to the Landau quantization, so that 
the magnetization and magnetoresistance oscillate by the 
frequency (/q.) of the area of the a orbit as a function of 
1/H. The high field enables electrons to tunnel the Bril- 
louin zone gap and to move along the larger closed orbit 
named as the /3 orbit. Thus, the oscillation frequency 
(//3) of the area of the P orbit gradually appears as H 
increases. This is the magnetic breakdown, which has 
been explained in a semiclassical theory by Falicpv and 
Stachoviakcl' based on Pippard's network model.EP 

In Falicov and Stachoviak theory, the oscillations with 
various frequencies (/„, //j, /2/3 = 2/^, ff^+a = ffS + fa, 
etc.) are allowed, since the closed orbital motion is pos- 
sible. However, the orbit named as the f3-a orbit in Fig. 
1 is prohibited, because the /3-a orbit is not allowed to 
execute the cyclic motion in the semiclassical picture. 
The oscillation corresponding to the area of prohibited 
orbits such as the (3-a orbit is called the quantum inter- 
ference oscillation. Shiba and FukuyamatP have shown 
the existence of the quantum interference oscillation in 



* E-mail: kishigi@sci.himeji-tech.ac.jp 



the transport in the case of the magnetic breakdown. It 
is now weJl-known as the Stark quantum interference os- 
cillation.LP On the other hand, it was believed that the 
quantum interference oscillation in the thermodynamics 
quantity dose not exist. 




Fig. 1. (a) The solid Hnes are the Fermi surface in our calcula- 
tions. The arrow indicates the direction of the electron orbital 
motion. For example, the effective masses for the a and /3 or- 
bits in (BEB|"R-TTF)2X are rria/mo '^ 1-4 and m^/mo ~ 3.8, 
respectively.B u^ (b) Their figures indicate the various orbits. 



ever, by using the simple tight-binding model 
3 one of the authors first found the quantum inter- 
ference oscillation, that is, the (3-a oscillation in the mag- 
netization M{H,n), whose frequency (//3-q — fp — fa) 
corresponds to the area of the (3-a orbit in Fig. 1 
by calculating the free energy, E{H,n), fuU-quantum- 
mechanically at the fixed total electron number, n. 

Recently, the existence of the (3-a oscillation in the 
magnetization was-sxperimentally observed in k-(BEDT- 
TTF)2Cu(NCS)2BI1J and a-(BKDT^TTE^sKHg(SCN)4.lli 

From other aspects of theoryJlil'EjEj'EZl'ta-' the exis- 
tence oLlhe (3-a oscillation in M{H,n) was^onfirmed. 
Nakanco' and Alexandrov and BratkovskyQ account 
for the (3-a oscillation in the magnetization M{H,n) 
by using two independent energy band model neglect- 
ing the magnetic breakdown. Harrison et. alli3) cal- 
culated by the density of states derived from Pippard 
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network ippdeLE' and found the (3-a oscillation. In addi- 
tion theytil'EjIla' reveal that there is no P-a oscillation 
in the magnetization M(H,fi) when the chemical poten- 
tial, /i is fixed. In other words, they clarify that the 
quantum interference oscillation such as the /3-a oscilla- 
tion is caused by the chemical potential oscillation. In 
fact in our model, in the case of the fixed chemical po- 
tential, the P-a oscillatiDn in M{H, /i) did not appear .E2P 
Moreover, Sandu et al.EZPand Han et alS3 found the P-a 
oscillation in M{H, n) from the full-quantum- mechanical 
calculation by use of more realistic tight-binding model 
based on the band calculation. n n n 

In our spinless tight-binding model,^^'^ the energy 
spectrum becomes so-called Hofstadter's butterfly dia- 
gramH in the presence of magnetic field. Various oscil- 
lations with /q, /^, /2^, //5-Q, etc., in M{H, n) are due to 
the chemical potential oscillating. In that calculation the 
oscillation with fp+a (P + a oscillation) is much smaller 
than P — a oscillation, although the P + a oscillation is 
obtained to be larger than, the P — a oscillation in the 
independent band model.M The origin of the suppres- 
sion of the P -\- a oscillation in the magnetic breakdown 
system is not clear yet, but since the P + a oscillation is 
large in the case of the fixed chemical potential,ll2P a can- 
celation between the effects of the magnetic breakdown 
and the chemical potential oscillation will cause the sup- 
pression of the P + a oscillation. If we take account of 
the spin, the energy spectrum will be divided into two 
energy bands. Even if the total number of electrons is 
fixed, the electron number for each spin is not fixed. As 
a result we may expect a richer effiect than the reduction 
factor in the semiclassical theoryEP by taking account of 
the effect of spin in the system with the magnetic break- 
down. However, it has never been studied that how the 
spin affects the "prohibited" quantum interference oscil- 
lation {P-a oscillation) and quantum magnetic oscilla- 
tions {a, P, P+a, 2/3 oscillations, etc.). Therefore, we 
need to investigate the effect of the spin on the magnetic 
breakdown and quantum interference oscillation. 

In the thermodynamics quantity such as the magneti- 
zation, P+a oscillation are observed in addition to a, 2a, 
/3-a, P and 2/3 osciUatians, in the quasi-two-dimensional 
organic conductors.EJ'H'F.T However, the peak at //j-q 
reported by Meyer et alS3 is very smalL-jwhile that peak 
is large in the experiment by Uji et alO Both de Haas 
van Alphen experiments measured at the same material, 
K-(BEDT-TTF)2Cu(NCS)2, are done at the same field 
region (22 T ~ 30 T) and the same temperature {^ 0.4 
K). The difference in the amplitudes of /3 — a oscilla- 
tion may be attributable to the effect of spin in a tilted 
magnetic field. 

In this paper, we calculate M{H,n) from the energy, 
E, at the canonical ensemble by adding the Zenmaii term 
to the previous spinless tight-binding modcl,&Bt2l' and 
try to analyze the effect of the spin on the magnetic 
breakdown and the quantum interference oscillation. 

§2. Formulation 

Let us study the two-dimensional tight-binding model. 
The lattice spacings in x- and y-axes are a and b' , respec- 



tively. The Hamiltonian at ff = is 

n^ic + v. (1) 

k.a 



*J>o' 



V = F^cosQ-r,C/^^a,,. 



(3) 



In the above a =t, |, CJ ^ is the creation operator of a 
spin electron at I site, e(k, cr) — — 2ta cos akx — 2tb cos b'ky 
and Q — (0,7r/6') is the lattice potential vector, where 
the Brillouin zone is —ir/a < kx < tt/o, — 7r/&' < ky < 
n/b' . In the momentum space, eq. (3) is written as 

^ ^ {CHkx,k'+{ '^^ 



V 



E 

.k'.rn.a 



1 



X C{kx,k'y + -—,G) + h.c.}, 



(4) 



where summation in k' and m are done in 
ky 



7r/26' < 
< 7r/26' and m = 0, 1, respectively. By V, the Bril- 
louin zone is halved so that the first Brillouin zone be- 
comes —ir/a < kx < 7r/a, — 7r/6 < ky < ir/b, where 
we choose b = 2b'. In this paper, we take parameters 
as 2/3-filling, 4/ia = 1-0 and V/ta = 0.09. Then, the 
Fermi surface is obtained as shown in Fig. 1, which has a 
characteristiolopology in quasi-two-dimensional organic 
conductorsE'cl' From this Fermi surface we can see that 
the possible cyclotron orbits are the closed a orbit and 
the breakdown /3 orbit. 

When the magnetic field is applied parallel to z axis, 
we take the Landau gauge A={Q,Hx,0), and we write K. 
by Peierls substitution as 



Hp— 



E 



{-2ta{cosa{k'^ 



nS)) - -gusHa} 



C^k'^+ n6, ky + m-,a)C{k'^ + nS, ky 
+ {-tbexp[ib'k'y + im'K]C^k'^ + nS, k'y - 



V 

■K 
'b' 



m-,c7) 



C (fc; + {n + 1)6, k'y + m-,a) + h.c.}, 



(5) 



where S = eaH/hc — (0/(/)o)(27r/a), cj) — ab'H is the flux 
passing through a unit cell and (po = 2T:hc/e is a unit 
flux, and h = (t>/4>o is the number of the flux quantum per 
unit cell. We represent magnetic fields by h, henceforth. 
When h is a, rational number, namely, h — 4>/ (pQ = p/q 
with p and q being mutually prime integers, the matrix 
size of Tip IS q y. q. In eq. (5), < fc^ < 5/p and 
Q <n < q. As a result, the matrix size of Tip-I-V becomes 
2qx2q. In Hp, /ie is the Bohr magneton and 5 ~ 2. If the 
Zeeman term in 7ip-|-V is excluded, this model becomes 
the previous spinless modelBoEEP 

Under the condition of the fixed total electron number, 
n, the ground state energy per site is calculated by 



1 



(6) 



where Ns is the total site number and Cj is one electron 
eigenvalue, and index j includes the spin index. The 
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authorSjBu'H and Harrison ct alllH) also have obtained 
the similar result. 



magnetization is given by M{h, n) = —dE{h, n)/dh. 

The magnetic field cannot be changed continuously in 
this formulation and the experimentally accessible field 
[H ^ 10 T) is difficult to study due to the large value of 
q. For exaixiple, as a = 12.71 and b = 8.4A for (BEDT- 
TTF)2xj3# h = p/q is about 1/400 when if ~ 10 T. 
We calculate at higher field {H - 100 T, /i - 1/40). 

In order to see the effect of the Zeeman term clearly, we 
define g as gh ^ i^gfi^H /ta- For example, in quasi-two- 
dimensional organic conductors with the transfer inte- 
gral, ta = 250A;B,Q'a^ and the lattice spacings mentioned 
above (a = 12. 7A and h = 8.4l) we get 5 ~ 5 ~ 2.0. 




10 20 30 40 50 60 1/h 70 



Fig. 2. Magnetization as a function of the inverse field. 



§3. Result and Discussion 

We show M{h,n) at 5 = 0,0.2,0.5 and 0.8 in Fig. 2. 
and the Fourier transform of M{h,n) in Fig. 3. 

First, we see Figs. 2(a) and 3(a) when g — 0. In 
Fig.3(a), these peaks (/« ~ 0.08, /2a ^ 0.16, Ua ^ 0.32, 
fp ~ 0.67, /2/3 ~ 1.33 and fp^a ^ 0.59) correspond to 
the area of a, 2a, ia, /3, 2/3 and /3-a orbitsc2P (the ratios 
of area of a and /? orbits to the first Brillouin zone are 
0.08 and 0.67, respectively). The oscillations of the long 
period (/„) and that of the short period (fp) in M{h,n) 
are seen in Fig. 2(a). As h increases, the amplitude of 
the oscillation with fp becomes gradually larger. This 
behavior is due to the magnetic breakdown. Although 
the f3-a oscillation in addition to a, (3 and 2/3 oscilla- 
tions exists, a peak for /3+a oscillation is very small. 
This property has been shown previously by one of the 
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Fig. 3. The Fourier transform amplitude (FTA) of M{h, n). The 
range of the Fourier transform is 10 < h^^ < 60. 



Next, from the Fourier transform of M{h,n) for g ^ 
(Fig. 3(b)~(d)), we can see that the peak at fp+a 
becomes large at g = 0.2, and this peak is the same 
order of magnitude as at fp, whereas the peak at fp-a 
is very small at g = 0.5. When g — 0.8, peaks at /„ 
and fp±a become very small, whereas peaks at /2a and 
//3±2a become large, where the (3 ~ 2a oscillation is the 
quantum interference oscillation like (3 — a oscillation. In 
order to understand the ^-dependence of these peaks, we 
show amplitudes of these peaks at /„, /2Q, fp, fp±a and 
fp±2a as a function of Ij from ^ = to 1.0 in Fig. 4. 

In semiclassical theory, the effect of the spin is repre- 
sented by the reduction factor ,EP Rs =^ cos{^pTigm/mQ), 
where p is for pth harmonics of the each frequency in 
M{h,n) and m (toq) are cyclotron effective mass (free 
electron mass). For example, when we choose m/mo = 2, 
p=l^ta = 250kB, a = 12. 71 and b = 8.4l, the Fourier 
transform amplitude (FTA) of M{h, n) is reduced by 
the factor |i?sp = |cos(7rg)|^. We show 7|i?sP in Fig. 
4, where we take 7 = 0.001. From Fig. 4, the g- 
dependences of the FTAs of /„, /2a, //3, fp-a and fp-2a 
decrease as 5 increases, and all of these except /„ increase 
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Fig. 4. The FTAs of some peaks {/„, /2a, fp, ffj±a and //3±2q). 
Thick Dotted Unes indicate 7|i?sP- 



upon further increasing 'g. These FTAs are in agreement 
with |i?sP quahtatively. It may be due to the different 
values of the effective masses and p that the position of 
the minimum of the FTAs of these frequencies are differ- 
ent. However, note that the g-dependences of the FTAs 
of fi3+a and fi3+2a are quite different from |i?sP; these 
ampUtudes first increase as g increases and decrease as 
g increases further. 

The anomalous ^-dependences can be understood as 
follows. In the spinless case (5 = 0), amplitudes oi P + a 
and (3 + 2a oscillations are suppressed by the cancelation 
of the magnetic breakdown and the chemical potential 
oscillation. By finite 5 the chemical potential oscillation 
is strongly suppressed. The suppression of the chemical 
potential oscillation is consistent with a strong reduction 
of /3 — a and j3 — 2a oscillations in small g as seen in Fig. 
4. On the other hand, the magnetic breakdown is not 
affected by g. As a result, (i + a and fi -\- 2a oscilla- 
tions becomes large as "g increases due to an incomplete 
cancelation of the chemical potential oscillation and the 
magnetic breakdown. 

When the field is tilted from the k^ axis to the k^- 
ky plane by 0, the component of the magnetic field for 
the Fermi surface along the kz axis is H cos 6. That for 
the Zeeman term is not changed by tilting of H. There- 
fore, the effect of magnetic fields for the Zeeman term is 
enhanced by tilting the field. 

As seen in Fig. 4 the amplitude of /3 — a oscillation 
depends on g, i.e., the tilting angle 6. Therefore, we ex- 
pect that the_different results between Meyer et al.EI 
and Uji et al.E3 are due to a direction of magnetic fields 
in both experiments. In fact, the frequency of ii oscil- 
lation in both experiments are not the sameOEj-' The 
relative magnitudes of the amplitude of each frequencies 
(/a, //3, //3-a, fp+a, etc.) are varied by changing g. 
Particularly, the oscillations with f/s+a and ff3+2a have 
the anomalous ^-dependences. Thus, we should consider 
the system of the magnetic breakdown with the effect of 
the spin, which h|a£_aot,beeii, iakeu into account in the 
previous studies.iBBtSBHy.y 



§4. Conclusion 

From the full-quantum mechanical calculation by us- 
ing the simple tight-binding model including an electron 
spin, we analyze the effect of the spin on the quantum 
interference oscillation such as P — a and /? — 2a oscil- 
lations and the quantum magnetic oscillation such as a, 
2a, (3, (3 + a and (3 + 2a oscillations. It is obtained that 
the relative magnitudes of all oscillations (a, f3, f3 + a, 
P — a, etc.) are changed by the Zeeman term. In par- 
ticular, the amplitudes oi p + a and P + 2a oscillations 
have anomalous gf^-dependences. 

We expect that the ^-dependences of the amplitudes 
of these oscillations in our result will be observed in the 
experiment of tilting magnetic field. 
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